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We consider two Jaynes-Cummings cavities coupled periodically with a photon hopping term. 
The semi-classical phase space is chaotic, with regions of stability over some ranges of the param- 
eters. The quantum case exhibits dynamic localization and dynamic tunneling between classically 
forbidden regions. We explore the correspondence between the classical and quantum phase space 
and propose an implementation in a circuit QED system. 
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I. INTRODUCTION 

The Jaynes-Cummings (JC) Hamiltonian is the canon- 
ical model for atom-light interactions, describing a single 
confined bosonic mode interacting with a two level sys- 
tem (qubit). This is sufficient to describe a wide range of 
phenomena in cavity Quantum Electrodynamics (QED). 
Systems of coupled JC cavities, the Jaynes-Cummings- 
Hubbard (JCH) systems , have been suggested for a 
diverse range of optical applications such as an optical 
analog for the Josephson junction [T] and Q-switching[2]. 
Networks of JC systems have also been predicted to ex- 
hibit phase transitions [3H5]. 

Improvements in the realization of photonic cavities 
in the lab have made possible exploration of Jaynes- 
Cummings systems 6-8] in the strong coupling regime 
in a variety of platforms. A current implementation of 
interest is in circuit QED, where a superconducting op- 
tical resonator is capacitively coupled to a Cooper-pair 
box. This is equivalent to a single cavity mode of the EM 
field coupling to a two level atom. The advantage of cir- 
cuit QED is that coherence times and atom-field coupling 
much greater than that can be achieved with visible and 
near infra-red systems. This makes circuit QED a po- 
tential medium for quantum computing, and already has 
been used to implement an 2 qubit Shor's algorithm[9]. 

The original proposals for observing quantum phase 
transitions in JCH systems [3HS] called for large numbers 
of identical systems. Constructing large arrays of cav- 
ities which are sufficiently coherent and identical poses 
a significant challenge. Exploiting long coherence times 
can allow some analogous effects to be studied by trading 
large-scale phenomena for small-scale, long time phenom- 
ena. For example, there is an isomorphism between the 
periodically kicked rotor and the Anderson tight binding 
model[10]. The Anderson model predicts localization for 
particles in a disordered lattice, and for dimension greater 
than three exhibits a second order phase transition be- 
tween metallic and super fluid phases. This has been 
recently demonstrated in the time-domain as a kicked 
system with cold atoms [TT]. 

We examine the dynamics of a pair of periodically cou- 
pled kicked JC systems using both quantum and semi- 
classical treatments. For two kicked coupled JC systems 




FIG. 1: (color online) Schematic of a possible supercon- 
ducting stripline cavity implementation of the kicked system. 
Transmon qubits are centered on each cavity at a and b with 
an atom-photon coupling ft, and with the inter-cavity cou- 
pling, k, controlled by an applied voltage at c. 



the semi-classical dynamics are non-integrable with a 
complicated phase space composed of regular and chaotic 
regions. The quantum case exhibits similar structure, 
which converges to the classical as the number of excita- 
tions in the system increases. 

Periodic systems, such as delta kicked rotors and tops, 
are widely used to study the link between classical and 
quantum chaos [T^]. Several interesting correspondences 
between the two regimes have been identified such as 
dynamic localization with regions of stability JS] and 
Lypanov exponents with entanglement generation[14j. 
There are many open questions about the nature of 
quantum systems with semi-classical dynamics that ex- 
hibit chaotic behavior, particularly in time varying 
systems [15]. 

We discuss a possible experimental implementation 
(figure [lj in a circuit QED system, compatible with the 
current state of the art and thus allowing an experimen- 
tal investigation of quantum chaos effects in a fast devel- 
oping field. Superconducting strip-line cavities coupled 
to transmons provide a JC couplingwell into the strong 
coupling regime [6 , and the architecture provides a sim- 
ple means for producing the kicked coupling (k) through 
an intermediate qubit [T^]. 
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II. MODEL 

The JC Hamiltonian, in the rotating wave approxima- 
tion is 

H JC = Ao- t t7 + / 9(crto + cro t ). (1) 

with <r(a) the atomic(bosonic) annihilation operator, A 
the atom-photon detuning, coupling energy (3 and we set 
h = 1. H JC commutes with the total excitation num- 
ber operator, L = a* a + crcr|17j. Therefore the total 
excitations in the cavity, I, is a good quantum number. 
In the bare basis, the eigenstates are 

= sm9i\g,l) + cos6»;|e, I - 1) , , 

|-,/) = cos 9i\g, I) -sin6i\e, I -1), yZ) 

where 

tan0z = 2M/(A + 2 X i), (3) 
H JC \±l) = (± X (l)-A/2)\±l) (4) 

and 

X (l) = V/3 2 / + AV4 

is the generalized Rabi frequency. Note the \fl depen- 
dence in interaction strength. The an-harmonic energy 
spectrum is the source of much interesting behaviour: In 
JC cavities it leads to photonic blockade [l8l fT9]. provid- 
ing an effective photon-photon non-linearity. In the sys- 
tem under consideration, the incommensurate energies 
result in dynamic localization, as will be shown below. 
The hopping term, 

K = K(a\a 2 + a) 2 ai) , (5) 

describes an interaction between the two cavity modes 
which allows photons to move from one to the other with 
hopping rate k, for example, via evanescent coupling in 
photonic crystals, or, in the case of circuit QED, capaci- 
tive or inductive coupling [20 . In our model the coupling 
is turned on periodically at times t = nT for a short 
duration r. Here, T is the period between kicks and n 
an integer. If r is sufficiently short (r <C 1//3), then the 
interaction can be described by a delta function "kick" : 

H = H( c + H 2 1C + S T K' (6) 

where Hf are the JC Hamiltonians for cavities 1 and 2, 
St is a periodic delta function with period T and K' — 
Kt. We also require that t ^> 1/w, so that the rotating 
wave approximation is valid. 

The three dimensionless parameters, kt, Tj3 and A/3, 
are sufficient to specify the dynamics of H . For simplicity 
we consider only the quasi-resonant case, A ~ 0, where 
the key features of the system are most easily elucidated. 
This makes sin^ = cos6*f = -h= in equation 2. 



The coupling term breaks the individual excitation 
conservation of each JC system, but commutes with the 
total L = Li + L 2l thus we can consider cases of total 
excitation number individually. For a single excitation, 
L = 1, the excitation oscillates between cavities trivially, 
with frequency kt, and so we do not dwell on this case. 
For all L > 1 we find rich behavior with signatures of 
quantum chaos. However, here we confine ourselves to 
L = 2 in the quantum case, and the semi-classical equiv- 
alent. Although the dimension of Hilbert space is just 
8, many of the features of quantum chaos are already 
present, and it is this case which will be most accessible 
experimentally. 

A. Semi-Classical dynamics 

We derive the classical equations of motion by taking 
the expectation value of the Heisenberg equations of mo- 
tion (see, for example, [2T] ). Between kicks each system 
evolves separately as 

(d) = E = -ifiS, 

(a) = S = iAS + if3ES z , (7) 
(d z ) =S Z =2fii(SE* -S*E), 

where E, the E-field, and S, vectors on the Bloch sphere 
are now classical quantities. For no detuning the uncou- 
pled equations of motion are equivalent to that of a pen- 
dulum with the momentum E and S z — cos 9, the height 
of the bob. This motion has two constants of motion, 

Ni = \E h2 \ 2 + i(S zi + 1) 

S*+4S*S = 1. {J 

While this has an analytical solution in terms of elliptical 
functions, in practice it is easier to numerically integrate. 
The kick is given by the map 

/ E t \ /cos k' sin k' \ ( E x \ , Q , 

{E 2 ) n+1 -[- S m K > co S K'){E 2 ) n - W 

The kicked hopping leads to non-integrable dynamics, 
so that the only constant of motion is now Ni + N 2 = N. 
In general this results in a chaotic phase space, however, 
for some values of k and T there will be regions in which 
the motion is semi-regular. These regions are described 
by KAM (Kolmogorov- Arnold- Moser) theory [25]. In an 
unperturbed system the path in the d dimensional phase 
space in action-angle variables lies on the surface of a 
d-torus. If the periods in each dimension are sufficiently 
incommensurate then the system is confined near a de- 
formed torus for small perturbations. The system be- 
comes increasingly chaotic as the perturbation is turned 
up, leading to destruction of some tori. The phase space 
is then a chaotic sea with islands of stability which are 
topologically separate, from the chaos as well as each 
other. Eventually the perturbation destroys all these re- 
gions and the dynamics are fully chaotic. 
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FIG. 2: Classical phase space strobe plots of the E-field in 
each cavity over 200 kicks from several initial points, a) kt = 
1.3, /3T = 0.1. In the small T limit the total total energy in 
the electric field, E\-\-E\ is stable, leading to non overlapping 
rings, b) kt = 0.4, f5T = 1.7. Phase space is mostly chaotic 
except for the four regions where the energy in the system 
is confined mostly to one cavity. As there are 4 degrees of 
freedom, and only a single constant of motion, plots of the 
electric field in each cavity do not convey the entire dynamics. 



The centers of stability that survive the longest are 
usually found around short periodic orbits. In this kicked 
system, however, there are in general no single-period or- 
bits, making the motion difficult to determine the precise 
point at which the phase space becomes fully chaotic. 
However, numerical simulations for the N = 2 case indi- 
cate that for small kt the most persistent KAM tori are 
around = V% sin (kt) 2 , N 2 ,i — v / 2cos(kt) 2 (Figure 
[2Jd). That is, in these four regions of phase space the 
energy in the system remains localized to a single cav- 
ity. Each period As kt is increased these regions become 
leaky (cantori) and eventually disappear, after which the 
phase space is fully chaotic. 

The value of kt at which the system becomes chaotic 
is dependent on T. The period for a small electric field in 
a cavity is 2tt; when f3T is resonant with this the KAM 
tori are destroyed with much smaller kt. Unlike other 
kicked systems, this system is still regular for some kt at 
the resonances due to the non-linear nature of the pertur- 
bation that each cavity sees. The range of parameters in 
which this mode occurs is shown in (figure [3^,) where the 
destabilizing effect of the resonances can be seen around 
(3T = 2nn. 

We can also consider the limit in which kt is larger 
then the kick period, f3T. In this limit the electric field 
decouples from the atomic degrees of freedom and the 
energy in the electric field oscillates between the two cav- 
ities (figure [2Jd) and we have separate regions which con- 
serve the total energy of the field. For small kick period, 
T <C (5 there is a center of stability around S z \ = S z2 = 0, 
dynamically confining the atoms to their ground states. 



B. Quantum Dynamics 

We find that the quantum dynamics exhibit some qual- 
itatively similar behavior to the classical case, however, 
there are also effects which arise which are specifically 
quantum in nature. 

To explore these dynamics we define the Floquet op- 
erator U f which evolves the system from time t — nT + 
to t = (n+ 1)T+: 

Uf = e ~KHi C +H^)T el K = e -iH Q T e iK (1Q) 

The dynamics of a kicked system can be studied though 
the eigenstates, fi of U. On application of U the Flo- 
quet states pick up eigenphase e lXi . Thus the problem is 
equivalent to a time invariant Hamiltonian. This allows 
the calculation of the long term behavior of the system. 

The quantum equivalent of KAM tori can be under- 
stood as dynamic localization|23j: States which are ini- 
tially in the localized regions have exponentially sup- 
pressed diffusion into chaotic areas of phase space. 

If some state tp is well represented by a small number of 
basis states, tp® we may consider i\> to be localized to some 
degree. This can be quantified with the participation 
number(P)[21]: 

P(i>)= (^El«#°)l 4 ^ _1 ( u ) 

which we have normalized by the total dimension d of 
the space. P is 1/d when 1(^1^)1 = 1 for some i and 1 
when ip projects evenly onto the ipf. One can consider 
this to be a indication of quantum ergodicitv|25). 

While P is dependent on the choice of basis (ie. we can 
always choose some basis with tp as a base), comparing 
the eigenstates of the unperturbed Hamiltonian to the 
perturbed best represents the degree of mixing [26], We 
therefore take the k = eigenstates as the basis, and 
increasing k leads to Floquet states with increasing P. 

Figure [3]d shows the average participation number of 
the Floquet states over a range of kt and f3T for a system 
with two excitations. We denote the subspace of states 
with two excitations in the one cavity as |-0?}s, and like- 
wise the states with one excitation in each cavity as tpjs. 
The regions where P is small corresponds to states with 
both excitations in the same cavity being dynamically 
separated from states with excitations in both cavities, 
ie. an approximate symmetry of Uf. 

The suppression is destroyed by resonances which oc- 
cur at T = i = ^,n(l + ^),n{l - ^), which are 
solutions to 

V2T = mT,m € L 

At these valuies the phase accrued after each period is 0, 
and so there is no destructive interference. This implies 
that is indeed dynamical localization suppressing disper- 
sion in the system. For example, when T — IL j^,the 
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FIG. 3: (color online) a) Classical: Average N2 over 1000 
kicks for points initially at E2 = v2- The red regions rep- 
resent parameters where most of the energy in the system 
is localized to a single cavity, ie. The chaotic border. The 
KAM tori are destroyed with relatively small kick strength at 
f3T — 2mr. b) Quantum: Average participation number of 
Floquet states. The regions of localization are qualitatively 
similar to those in the classical case, however, note that the 
T scale is different in each graph. Vertical lines mark the 
location of resonances. 



states in \ip 2 ) pick up no relative phase to states with 
E = 0. This removes the interference suppressing trans- 
mission into these states, and destroys the localization. 

In figure [3Jd) we can see, for the atomic limit, that the 
dependence of localization on the parameters correspond 
qualitatively to the semi-classical case, though with im- 
portant differences. The frequency at which the classical 
cavities oscillate depends continuously on the energy in 
the cavity, and in general is different from the Rabi fre- 
quency of the quantum case; these two only coincide in 
the limit I —¥ 00. Thus, the locations of resonances are 
different in the two regimes. 

Note also that in contrast to the classical case, the 
resonance removes the localization for arbitrarily small 
kt. Resonances in the classical case are not sharp, due 
to the energy dependent frequencies. 

For time independent systems, chaos can be studied 
via the statistics of energy levels, however, in periodic 
systems, the eigenphases of the unitary operator are not 
observable. Ergodicity of can be explored experientially 
by comparing the expectation of observables in the sys- 
tem to an ensemble of random states. For a chaotic sys- 
tem, the unitary map Uf has no symmetries, and so we 
expect the average state to be no different from a random 
one chosen with the appropriate measure. Figure [4] shows 
the long-time average of some experientially observable 
quantities, and the expected average of a random state. 

Classically, islands of stability are topologically sepa- 
rated, forbidding transitions between them. Quantum 
dynamics admit such flow of probability in phase space 
by a mechanism called dynamic tunneling and has been 
observed experimentally in a variety of systems [27]. Al- 
though this mechanism is distinct from the usual tun- 
neling, as there is no potential barrier to overcome, the 
system nevertheless moves across classically forbidden re- 



FIG. 4: Average participation number and long term aver 
ages of measurable quantities (mean expectation of an en 
semble of random states) over increasing kt with f3T 
1.2. Circles(Dashed): (a^ zl + l)/2, Diamonds(Dots 
Squares(Dot-Dashed):<r z icr Z 2, Triangles: Average participa- 
tion number(Dot-Dashed) . 
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FIG. 5: (color online) Evolution of the system initially in \ip 2 } 
over 1000 kicks with f3T = 1.2 and a) kt = 0.1 b) kt = 0.2. 
Blue line is the expectation of excitations in cavity 1. Purple 
line is the expectation of finding both excitations together in 
the one cavity. 



gions in phase space. 

In the k = limit there is a two fold degeneracy for all 
Floquet states due to the H{ c , if/ symmetry. A state, 
initially in l^j) in cavity one is in a superposition of 
two Floquet states, | ± / 2 ) , which have equal projections 
in both cavities, but still in the \ipf) subspace: 



1 



(l + / 2 ) + |-/ 2 » 



(12) 



The perturbation breaks the degeneracy, leading to an 
approximate separation in the eigenphases, 4>. Each kick, 
the two Floquet states composing \ip) are separated by a 
phase-angle of </>. After ^ kicks the phase separation is n, 

and \tp) has evolved to the state -^(|+ / 2 ) -|- f 2 ) = l^f), 

i.e. completely in the other cavity. Figures [5^,) andjs]}) 
show the transmission between the two separated local- 
ized states for kt = 0.1 and kt — 0.2 respectively. The 
two excitations in the system oscillate between cavities, 
though are strongly localized to the ijji subspace. As 
kt increases so does <fi, and the localization to the \ipf) 
subspace decreases. 
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III. EXPERIMENTAL IMPLEMENTATION 

While the effects discussed apply to any implemen- 
tation of JC systems, circuit QED (cQED) presents it- 
self as one of the most viable platforms due to the large 
coupling coefficients and long coherence time, relative to 
other cavity QED systems. 

Current experiments in cavity QED, where a transmon 
is coupled to a resonating microwave cavity, have charac- 
teristics which could allow a successful realization of this 
kicked system. A cQED setup with lj/2it = 6.92GHz, 
P/2tt = 347MHz and coherence time of order 1/is has 
been achieved recent ly [HI H] • 

The localization transition occurs around act « .1 and 
for the delta-function kick approximation to be valid we 
need the pulse time t <C 1/(3. For the coupling strengths 
cited above, this requires a pulse time of t m lCT 10 s and, 
therefore, k order 1GHz. Between pulses K must be of the 
same order as the decoherence rate (ie. ~ 1MHz) such 
that the dispersion due to the constant inter-cavity cou- 
pling is small over the time of the experiment. Thus a se- 
quence of ~ 100 kicks could be applied within the coher- 
ence time. We have seen that this is long enough to ob- 
serve dynamic tunneling and localization/dclocalization 
by inlcuding the decoherence and dephasing explicitly in 
the simulation. 

The tunable hopping term could be achieved using an 
intermediate qubit coupling such as in [THl HO] • In such 
scheme's the effective coupling is of order 

where /3i3, P23 and A 3 are the coupling strengths of each 
resonator to the intermediate qubit and it's detuning re- 
spectively and A3 3> /3. This requires the coupling to 
the intermediate qubit to be significantly greater than 
the other couplings. The detuning can be controlled in 
situ, allowing the coupling to be switched on and off. 

Spectroscopic measurements can be used to determine 
the final state 0- Although there will be significant in- 
teraction with the environment, the only final states of 
interest are those that still have two excitations. One 
can therefore largely remove the effects of atomic re- 



laxation and photon dissipation with a post-selection 
scheme, given a temperature smaller then then the char- 
acteristic energies of the system. De-phasing terms will 
still be relevant, however, these are generally ignorable 
over the time frames considered [5] • 



IV. DISCUSSION/CONCLUSIONS 

The phenomena discussed have been observed in other 
systems, such as dynamic tunneling and localization in 
cold atoms [27, 28\. Circuit QED allows direct control 
over many system parameters and direct measurement of 
the state of the system. This can be used, for example, 
to study the effect of noise by controlling the detuning 
parameter in situ. 

As circuit QED is proving to be an important field, 
with a wide range of possible applications, understand- 
ing chaotic behavior in these systems will be crucial. An 
experimental realization of the system seems quite pos- 
sible, although it is not without challenges, specifically 
in achieving a sufficiently large inter-cavity coupling. It 
would allow the study of the rich behavior that can be ex- 
pected in coupled Jaynes-Cummings systems, and open 
up new regimes for investigating quantum chaos. 

We have presented a simple model which exhibits a 
transition from localization to ergodicity and dynamic 
tunneling. Importantly, we see this behavior even for 
small Hilbert space dimension, which, although interest- 
ing behavior can be seen for any number of excitations 
above two, the lowest case most clearly conveys the as- 
pects we have emphasized. Furthermore, the two exci- 
tation case will most likely be the easiest to implement 
experimentally Constantly improving control in circuit 
QED systems means that it will be possible to study the 
higher dimensional cases. This could potentially allow a 
novel means for probing the transition between classical 
and quantum chaos. 
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